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ABSTRACT 

Context. The DANCe survey provides photometric and astrometric (position and proper motion) measurements for approximately 2 
million unique sources in a region encompassing ~80 deg^ centered on the Pleiades cluster. 

Aims. We aim at deriving a complete census of the Pleiades and measure the mass and luminosity functions of the cluster. 

Methods. Using the probabilistic selection method described in |Sarro et al.|pOT4) , we identified high probability members in the 
DANCe (i >14 mag) and Tycho-2 (V <12 mag) catalogues and studied the properties of the cluster over the corresponding luminosity 
range. 

Results. We find a total of 2 109 high-probability members, of which 812 are new, making it the most extensive and complete census 
of the cluster to date. The luminosity and mass functions of the cluster are computed from the most massive members down to 
~0.025 Mq. The size, sensitivity, and quality of the sample result in the most precise luminosity and mass functions observed to date 
for a cluster. 

Conclusions. Our census supersedes previous studies of the Pleiades cluster populations, in terms of both sensitivity and accuracy. 
Key words. Proper motions, Stars: luminosity function, mass function. Galaxy: open clusters and associations: individual: M45 


1. Introduction 

As one of the closest and youngest open clusters, the Pleiades 
has been extensively studied over the past century. From stellar 
activity to multiplicity and lithium abundance, many topics have 
been addressed using this fascinating cluster as a cornerstone. 
One very important subject is the shape of the origin of the initial 
mass function, a key tool in understanding the formation and 
evolution of the stellar populations in different environments. Its 
computation requires an accurate description of the stellar (and 
substellar) census. 

As a matter of fact, the study of the Pleiades as a group 
started much earlier. Both Claudius Ptolemy and Tycho Brahe 
identified a few dozen stars. Even Galileo Galilei did not do 
much better when he made the first telescopic observation of the 
cluster. Coming back to the 20th century, the stellar census ex¬ 
panded rapidly due to the increase in sky coverage and the depth. 
To name a few that are important because the naming conven- 

* Based on service observations made with the William Herschel 
Telescope operated on the island of La Palma by the Isaac Newton 
Group in the Spanish Observatorio del Roque de los Muchachos of the 
Instituto de Astrofisica de Canarias. 


tion and prefixes: |Trumpler| (|1921a]l and |Trumpler| (|1921b]l for 
Tr;|Hertzsprung|(|1947| l for ' HII ";|Artyukhina|(|1969 for "AK"; 


|Haro et'aL 


for "HCG"; [van Leeuwen et ah 


([T986 


___ i_Ll_. 

"Pels";|Stauffer et al. 

HHJ; [Pinfield et al.|(|2000irfor "BPL"; and [Deacon & Hambly 


P004 1 for "DH". The Pleiades has also been a privileged target 


for the search and study of substellar objects. The first confirmed 
brown dwarfs were in fact discovered in the Pleiades ( |Rebolo| 
|et al.||19^ |Basri et al.|1996) l, extending the mass function be¬ 
yond the hydrogen-burning limit. Several studies subsequently 
identified a growing number of substellar Pleiades members, 
providing hints that the mass function flattens or even increases 
at the very low mass end (e.g., Moraux et al.|2003 Lodieu et al. 


2012| and references therein) and suggesting different (or addi¬ 


tional) formation mechanisms. 

[Stauffer et al.| ( |2007| l have recently presented an exhaustive 
list of candidate and confirmed members that come from more 
than a dozen independent surveys of the Pleiades down to the 
limit of sensitivity of the 2MASS 6x Atlas and Catalogs. The 
total number of members and candidate members reported in 


their work adds up to 1 471 objects. More recently, Lodieu et al. 


\2Q\2) have used the UKIDSS GCS survey to search for low 
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mass members. The limited astrometric accuracy of the UKIDSS 
survey and the inconsistent probability membership calculations 
used in their study resulted in an incomplete and contaminated 
list of members ( Sarro et al.|2014| l. 

In |Sarro et al. ( 2014| l, we presented a new probabilistic 
method based on a multivariate analysis to derive reliable mem¬ 
bership probabilities using the DANCe catalog of |Bouy et al.| 
( 2013| l. In this article, we present a new release of the DANCe- 
Pleiades catalog including new photometry and updated mem¬ 
bership probabilities (Appendix [A), as well as a new analysis 
of the Tycho-2 catalog using the same technique (decribed in 
Appendix |^, in order to also tackle the high mass end of the 
cluster. We used them to derive an improved list of high proba¬ 
bility cluster members, described in section]^ Since the Pleiades 
has played a very important role in the history of brown dwarf 
searches, we assess different surveys again by looking for sub- 
stellar members in Section]^ In section [6T| we derive the lumi¬ 
nosities and effective temperatures of the Pleiades members. In 
section]^ we compute the luminosity and mass functions of the 
cluster and compare them to theoretical predictions. Finally, in 
sectionj^we provide emprical isochrones of the Pleiades single¬ 
star locus in the Bj,VT,u,g,r,i,z,Y,J,H, and Ks filters. 


2. Building the sample of members 

Deriving a reliable luminosity function representative of the true 
cluster’s population requires a good understanding of the selec¬ 
tion biases. In an effort to cover the entire luminosity range, we 
combined candidate members from two complementary sources. 

- Middle and low luminosity end: the DANCe survey is ex¬ 
pected to be complete over the dynamic range between 
14< i <23 mag (see Fig. EB- A new (second) release 
of the DANCe-Pleiades catalog ( |Bouy et al.||2013| is pro¬ 
vided with the present article, and described in Appendix [A| 
The main improvement consists in the addition of gri pho¬ 
tometry from the AAVSO Photometric All-Sky Survey DR7 
(APASS), as well as our own new Y-band photometry for 
a small subsample of candidate members at the bright end, 
when no measurement was available. These two additions al¬ 
low us to recover many members previously missed because 
of the lack of a good training set with complete photometric 
coverage at the bright end of the DANCe survey. In |Sarro| 
|et al.| ( |20I4[ ), we discussed various selection biases in the 
DANCe-Pleiades catalog relative to the choice of the vari¬ 
ables used for the membership probability calculations. As 
recommended in this article and for the rest of the present 
study, we use the membership probability defined using the 
RF2, and a selection threshold at probability greater than or 
equal to 75%. A total of 2010 sources were selected. 

- High luminosity end: the Tycho-2 survey provides astromet¬ 
ric and photometric measurements for 2.5 million sources 
brighter than V <12 mag. We complemented the Tycho-2 
Bj, Vt photometry with APASS, 2MASS, and CMC-14 pho¬ 
tometry, and applied the method described in |Sarro et'af] 
( 2014| l to the corresponding catalog within the same area 
as the DANCe survey. The selection process - and in par¬ 
ticular the choice of variables and selection threshold - are 
described in more detail in Appendix [B| We identified 207 
candidate members, which overlap with the DANCe selec¬ 
tion mentioned above. 

The DANCe and Tycho-2 samples overlap over more than 
3 mag and include 110 common sources. The overlapping region 
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Fig. 1. Sample size (not including the Tycho-2 analysis) as a function 
of the membership probability threshold. The raw (cyan) and corrected- 
for contamination (red) values are indicated. A dotted line indicates the 
0.75 level chosen for the present study. 


is believed to be incomplete because it corresponds to the limit 
of sensitivity of the Tycho-2 catalog and to the saturation limit 
(largely inhomogeneous over the field covered) of the DANCe 
survey. Even though the final catalog seems to cover the entire 
luminosity range, the domain between 12< i <14 mag must 
suffer from biases and incompleteness that are hard to quan¬ 
tify. Finally, a visual inspection showed that Alcyone and Electra 
had been missed by our selection. Both are indeed absent in the 
Tycho-2 catalog. Eor completeness, we added them to the final 
list of members, which include a total of 2 109 high probability 
members. 


Eigure[T] shows the sample size - raw and corrected for con¬ 
tamination - as a function of the membership probability thresh¬ 
old for the DANCe sample alone. The contamination rate was 
estimated using simulations as described in |Sarro et al.| ( 20T4) l. 
It shows that the total number of members after correction for 
contamination slowly increases to reach ~2 000 - 2 100 members 
at a threshold of ~0.2 - 0.3, and then rapidly increases at lower 
membership probabilities. The simulations probably largely un¬ 
derestimate the contamination level for the lowest membership 
probability values, explaining the sharp increase. But Eig.[2nev- 
ertheless suggests that between 200 and 300 sources below the 
conservative 75% threshold used in our analysis could be gen¬ 
uine members whose probability was lowered because of prob¬ 
lematic photometric measurements or very large astrometric un¬ 
certainties. As such they deserve attention and follow-up obser¬ 
vations to confirm their membership in the cluster. The same ap¬ 
plies to the Tycho-2 sample, in which ~10 - 20 members might 
have been missed by the 0.48 membership probability selection 
threshold (see Appendix |^. With this limitation in mind and 
to minimize the contamination level, for the rest of the present 
analysis we use the recommended 0.48 and 0.75 thresholds for 
the Tycho-2 and DANCe samples, respectively. 
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Fig. 2. Venn diagram representing the sample of candidate members 
identified in the DANCe survey (blue), in lStaulfer et al.|p()07^ (green), 
and in |Lodieu et al.| ( [2012| l (red, probability greater than 50% in their 
analysis). In each case the total number is indicated in bold face outside 
the circles. 


3. Comparison with previous surveys 


We compare the result of our selection with the list of 1 471 can¬ 
didate members from Stauffer et al.| ( |2007| l and 757 candidate 
members from [Lodieu et al.|P012| l. These two studies cover a 
similar area, and their dynamic ranges largely overlap, resulting 
in a total of 644 candidates in common. The two lists combined 
add up to 1 584 unique candidate members, of which 1 297 have 
a counterpart within 2" in the DANCe catalog. Figure shows 
a Venn diagram comparing the three samples of candidate mem¬ 
bers. A total of 812 candidate members identihed in the present 
study are not included in [Stauffer et al. ( [2007 1 or Lodieu et al. 
( |2012| l. Figureshows a (/, i - K) color-magnitude diagram for 
the 1411 candidate members with both i and K photometry. 

We investigated the candidates from Stauffer et al. ( |2007| l 
and jLodieu et al.| ( [2012| l with a counterpart in the DANCe survey 
but a membership probability below our 0.75 threshold. We ex¬ 
plored the Malahanobis distances in the hepta-dimensional space 
of proper motions, luminosities, and colors. Figure ^presents the 
results. Four Malahanobis distances were explored; 


- MDl: minimum distance to any of the Gaussian components 
in the bidimensional proper motion space; 

- MD2: minimum distance to the principal curve (isochrone) 
of the single star locus in the penta-dimensional space made 
of colors and luminosities; 

- MD3: minimum distance to the principal curve (isochrone) 
of the equal-mass binary star locus in the penta-dimensional 
space made of colors and luminosities. 

Each MD therefore tells whether a source has been rejected 
mostly because of its inconsistent proper motion (MDl) or pho¬ 
tometry (MD2 and MD3) or both. Figurej^shows that in general 
and as one could expect, objects with low membership probabil¬ 
ities have large Malahanobis distances in several spaces, while 
high probability cluster members are generally located at short 
Malahanobis distances in all spaces. A source with an incon¬ 
sistent proper motion (large MDl) is indeed most likely unre¬ 
lated to the cluster so is likely to display colors and luminosities 
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Fig. 3. i vs i-K color-magnitude diagram of the DANCe survey. Sources 
with membership probabilities greater than or equal to 75% are indi¬ 
cated as large colored dots. The color scale represents the membership 
probability. 


that are inconsistent with the cluster single and binary sequences 
(MD2 and MD3). These results suggest that in general the rejec¬ 
tion is robust because it is consistent in several dimensions. 

A number of members from jStauffer et al.] ( |2007| ) have mem¬ 
bership probabilities below our 0.75 threshold but above 0.1. As 
discussed above, these could well be genuine members as well. 
Additionally, a number of candidate members from jStauffer et al.] 
( 2007| l fall in the domain of incompleteness between the Tycho- 
2 and DANCe analyses, are missed by the present study, and 
contribute to the difference illustrated in the Venn diagram. 


4. Brown dwarfs in the Pleiades 


The Pleiades cluster has been one of the most popular brown 
dwarf hunting grounds. The study by Lodieu et al. ( 2012| l is dis¬ 
cussed above. In the following we compare the list of candidate 
substellar members of the Pleiades from Moraux et al.|(|2003]l. 


Bihain et al. ( 2006| ), Casewell et al. (2007 i, and Zapatero Osorio 
et al.j (|2014b i to the DANCe list of high probability members. 


Earlier studies exist but are superseded by the latest ones. 


4.1. \Moraux et al.\(2003l 

Based on proper motion (/, Z) photometric measurements, 
[Moraux et al.j ( [2003[ l have discovered 37 brown dwarf candi¬ 
date members of the Pleiades. The entire dataset of images used 
by these authors was included in the DANCe survey, ensuring a 
complete spatial overlap with the corresponding study. Only 31 
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4.2. \Bihain et af] \2006^ 


4.3. Casewell et al. 2007) 


Casewell et al.| ( [2007| l combined the UKIDSS ZYJHK images 


with the far red I, Z CFH12K images presented in |Moraux et al.| 
Article number, page 4 of|20| 


( |2003| l to search for new members. A total of 16 out of 23 candi¬ 
date members have a counterpart in the DANCe catalog. Five out 
of these 16 candidate members have very low membership prob¬ 
abilities. All but one have proper motions that are inconsistent 
with the cluster’s motion. The last one (PLZJ-56) has a proper 
motion consistent with the cluster’s mean motion but optical lu¬ 
minosities and colors inconsistent with the cluster’s sequence: 
the source appears to be slightly bluer (by ~0.2 mag) than the 
Pleiades sequence in all color-magnitude diagrams. This object 
could be a genuine member with a peculiar color, so it deserves 
further attention. Follow-up observations are needed to conhrm 
its membership status. The remaining 11 have high membership 
probabilities(>99%) and are considered to be good candidate 
mem bers. Table [T|giv es a summary of the 16 candidate members 
from [Casewell et ar] ( 2007] l with a counterpart in the DANCe 
survey. 


4.4. jZapafero Osorio et al] (2014b^ 


Fig- 4 . Distributions of Malahanobis distances for the sample of jStauff^ 
|et al.| ( [2007) l. The color scale represents the membership probability as 
reported in the DANCe catalog. 


candidates have a counterpart in the DANCe catalog within 2". 
A visual inspection of the CFH12K, MegaCam, and UKIDSS 
images shows that three of the six remaining (CFHT-PLIZ-5,14, 
18) are not detected in the CFFI12K, MegaCam, or UKIDSS im¬ 
ages, so we interpret them as false detections. The other three 
(CFHT-PLIZ-23, 34, 37) are fast-moving sources that visibly 
move by several arcsec over the 13-year period and can be re¬ 
jected as well. 

A total of 13 objects classihed as members by jMoraux et al.j 
( [2003 | l have membership probabilities below 10% and are con¬ 
sidered as contaminants. All of them have such very low mem¬ 
bership probabilities mostly because of proper motions inconsis¬ 
tent with the cluster’s mean mo tion. Table [T]gi ves a summary of 
the 19 candidate members from Moraux et al.| ( 2003] l rejected in 
our analysis. 


jZapatero Osorio et al.j ( |2014b) l have recently reported the dis¬ 
covery of 44 ultracool dwarf candidate members of the Pleiades 
based on proper motions and near-infrared photometry. Thirty 
two sources from their Tables 1 and 2 have a counterpart in the 
DANCe catalog ( Sarro et ar]|2014 and Table [^. Twenty have 
a membership probability close to 0 because their proper mo¬ 
tions are inconsistent with the cluster’s mean motion within the 
uncertainties. 

In most cases, the DANCe dataset used to derive these proper 
motions are excellent, with a long time baseline, a relatively 
large number of epochs, and reduced-;^^^ close to unity. The deep 
multi-wavelength photometry included in the DANCe catalog is 
also useful for testing their membership. Several are detected in 
either u, g, or r or in the old DSS plates, which is highly un¬ 
likely for such ultracool Pleiades brown dwarfs. All are detected 
in either the r or i-band with luminosities inconsistent with the 
empirical Pleiades sequence, as shown in Fig. On the basis of 
our more precise proper motion measurements, optical photom¬ 
etry, and probabilistic membership analysis, we rejected these 20 
candidates and classihed them as contaminants. Only 12 candi¬ 
dates of the 32 with a counterpart in the DANCe catalog have 
high membership probabilities according to Sarro et ari(2014|l. 


[Bihain et al. (2006|l complemented the CFH12K data presented 


in Moraux et al. ( 2003| l with deep near-infrared J-band images to 
search for new substellar members. A total of 26 out of 34 very 
low mass and substellar members of their list have a counterpart 
in the DANCe catalog. None of the eight that are missing are de¬ 
tected in our deeper MegaCam stacked images even though they 
are supposed to be brighter than the limit of detection, so we 
consider them as false detections. Finally, four of the twenty-six 
candidates with a counterpart in the DANCe catalog have mem¬ 
bership probabilities close to zero essentially because of incon¬ 
sistent proper motions, and are also considered as contaminants. 
The remaining 22 have high membership probabilities(>99%) 
and are considered as good candidate mem bers. Table [T|g ives a 
summary of the 12 candidate members from jBihain et al.| ( |2006| l 
rejected in our analysis. 


All 12 had been previously discovered and identified in the liter¬ 
ature or in the DANCe survey. 

Finally, half of the 12 sources without counterparts in jSarroj 
et al. (2014|l have a counterpart in the complete DANCe cata¬ 


log hut no reliable proper motion, so they were not included in 
the [Sarro et ah] ( [2014| l analysis and catalog. The DANCe com¬ 
plete catalog includes photometric information that helps dis¬ 
card a few more candidates. Calar Pleiades 31 is detected in the 
r band with 22.66+0.12 mag - too bright for a J=21.37 mag 
Pleiades ultracool dwarf. Calar Pleiades 23, 28, 29 and R2, R3 
are detected in the /-band with luminosities inconsistent with the 
empirical Pleiades sequence (see Fig.|^. Calar Pleiades 22, 24, 
25, 26, 30, and 31 do not have any counterpart in the DANCe 
complete catalog. According to the authors, Calar Pleiades 26 
is associated to a faint source in the UKIDSS Ks image, which 
they use to derive a proper motion measurement. A visual in¬ 
spection of the UKIDSS and WISE images shows that the cross¬ 
identification is wrong, because the counterpart is located 174 
north of the reported position of Calar Pleiades 26 and is asso¬ 
ciated to an optical source (dance-j03451605+2335088, within 
072). The corresponding astrometry and photometry is therefore 
wrong. The detection seen near (but not exactly at) the location 
of Calar Pleiades 26 in the same Omega-Prime images as the one 
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Fig. 5. i vs i - J color-magnitude diagram of the DANCe survey. High 
probability (>0.95) sources are indicated as cyan dots. Candidate mem¬ 
bers from Zapatero Osorio et al.| ( |2014b^ with a i-band counterpart in 
the DANCe survey are represented with red squares. Their abbreviated 
names are also indicated. Candidate members fro m [Zapatero Osorio| 
jet al.j ( |2014b^ previously discovered or reported in jSarro et al.| ( |2014] l 
are overplotted as large blue circles. 


used by these authors is extremely faint and, if real, is blended 
with the DANCe source 5148179. For the rest of the study, we 
consider Calar Pleiades 26 as a contaminant. A total of five can¬ 
didates remain, leading to a contamination rate of >61% in the 
Zapatero Osorio et al. p014b|l survey. 


Zapatero Osorio et al. 


From this original sample, 

( |2014a| l have recently obtained near-infrared spectra of 
the candidate members dance-j0347467B-i-2535166, dance- 
j0345333'3-b2334342, dance-j03481565H-2550089, dance- 

j0345584B-i-2341539, and Calar Pleiades 21 and 22. The 
spectroscopic features support their cluster membership. They 
also obtained an improved upper limit on the Z-band photometry 
of Calar Pleiades 25, which is consistent with membership but 
does not allow drawing a hrm conclusion about its membership. 


4.5. Importance of inspecting the images 


Proper motion surveys based on a couple of epochs have been 
commonly used to search for members of young clusters. The 
present analysis hnds contamination rates of at least 50%, 35%, 


43%, and 61% in Moraux et al. P003 i, Bihain et 

al.[ ([2006|l, 

Casewell et al.|(2007|l, and 

Zapatero Osorio et al. ( 

2014b 

1, re- 


spectively. The vast majority of identified contaminants have dis¬ 
crepant proper motions according to the DANCe measurements. 


which illustrates the need to have more than just two epochs to 
derive robust proper motions from ground-based images. 

Additionally, a significant fraction of the contamination is 
related to false detections and illustrates the need to always vi¬ 
sually inspect the images and not rely just on the catalogs. The 
latter inevitably contain errors and false or problematic detec¬ 
tions. These recommendations also apply to the DANCe cata¬ 
log itself, even though it is based on a larger number of images 
and epochs (so is generally more reliable). For that reason, we 
make the MegaCam r,i and UKIDSS Ks mosaics used for the 
DANCe survey available on the internet at the following address: 
http://visiomatic.iap.fr/pleiades/ The images can be 
browsed in an intuitive, interactive way thanks to the VisiOmatic 
tool ( jBertin et al.|2014| l. The Pleiades members identified in the 
DANCe survey are indicated, and their photometry, spectral en¬ 
ergy distribution (SED), best-fit effective temperatures (see Sec¬ 
tion |6.1| l, proper motions, and membership probabilities are re¬ 
ported. 


5. Empirical Pleiades isochrones 

5. 1. Empirical isochrones 

The robust membership probabilities and multiwavelength pho¬ 
tometry covering the spectral range from the near-UV the near- 
infrared offers a unique opportunity to measure precise empirical 
isochrones. Such isochrones prove extremely useful as a refer¬ 
ence for the study of other groups and clusters and can be used 
as an absolute empirical calibration for the evolutionary models. 
Table gives the apparent magnitudes for the Pleiades empiri¬ 
cal sequence in the u,g,r,i,Y,J,H,Ks over the dynamic range of 
the DANCe dataset. The f-band was used as reference, and the 
other luminosities were imputed by fitting a principal curve to 
the ridge of the empirical sequence in the corresponding color- 
magnitude diagram and by applying a manual offset (typically 
between 0.1 and 0.2 mag in color) to match the lower edge, as 
illustrated in Figure]^ To minimize systematics and errors, only 
sources with a membership probability greater than or equal to 
0.99 were used. Similarly, Table[^gives the apparent magnitudes 
for the Pleiades empirical sequence in the Bj, Vt, J, H, and K- 
bands obtained from the reanalysis of the Tycho-2 data and us¬ 
ing Vt as reference. Increasing photometric errors towards the 
faint end result in larger uncertainties for the location of the se¬ 
quence. The isochrones should therefore be considered and used 
with caution at the faint end (see Fig[^. 


5.2. Comparison with evolutionary models 

Figure |7] compares the observed sequence of Pleiades members 
in various color-magnitude diagrams to the models of [Allard] 
P014) l and |Bressan et al. ( 2012| l for an age of 120 Myr. Members 
with a probability greater or equal to 0.99 are represented to min¬ 
imize the contamination and associated noise in the comparison. 
The match is in general relatively poor in any color-magnitude 
diagram made of optical and near-infrared luminosities. As al¬ 
ready known (see, e.g., [Bell et aT]|2012[ [Baraffe et al.||199'8| 
and references therein), the difference between observations and 
models is larger at the faint end and suggests that the current de¬ 
scription of the atmospheres of late-M and L-type brown dwarfs 
presents a number of problems. Surface inhomogeneity (active 
zones, spots, Stauffer et al.|2003) Kamai et al.|2014 1 , and rota¬ 
tion ([Somers & Pinsonneault|2014 Kamai et al.|2014|l, molec- 


ular opacity (e.g., Faherty et al.|2012 Manjavacas et al.|2014| l, 

weather ( [Crossfield et al. 2014| l, and abundances (specifically the 
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Fig. 6. Empirical isochrones in the u,g,r,i,Y,J,H,Ks using the i-band as reference (colored curves) overplotted on top of high probability members 
(proba>0.99). 
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Table 2. Empirical isochrones of the Pleiades. More detailed version 
available on-line 


i 

u 

g 

r 

z 

Y 

J 

H 

Ks 

9.4 



9.36 



8.779 

8.617 

8.556 

9.9 


10.33 

9.886 



9.195 

8.962 

8.886 

10.4 


10.95 

10.42 



9.604 

9.301 

9.209 

10.9 


11.59 

10.97 



10.0 

9.629 

9.522 

11.4 


12.27 

11.53 



10.39 

9.942 

9.821 

11.9 

14.78 

13.0 

12.12 



10.75 

10.25 

10.11 

12.4 

15.77 

13.76 

12.74 



11.11 

10.55 

10.39 

12.9 

16.69 

14.52 

13.38 



11.45 

10.86 

10.67 

13.4 

17.52 

15.26 

14.04 

12.98 


11.79 

11.18 

10.97 

13.9 

18.25 

15.97 

14.72 

13.36 


12.12 

11.52 

11.28 

14.4 

18.93 

16.66 

15.39 

13.76 


12.46 

11.88 

11.62 

14.9 

19.59 

17.32 

16.05 

14.18 

13.32 

12.82 

12.26 

11.98 

15.4 

20.23 

17.96 

16.69 

14.62 

13.71 

13.2 

12.66 

12.36 

15.9 

20.87 

18.58 

17.3 

15.07 

14.12 

13.6 

13.07 

12.76 

16.4 

21.55 

19.2 

17.91 

15.53 

14.55 

14.0 

13.48 

13.15 

16.9 

22.24 

19.85 

18.52 

15.99 

14.98 

14.39 

13.88 

13.54 

17.4 


20.52 

19.16 

16.44 

15.37 

14.76 

14.25 

13.89 

17.9 


21.23 

19.84 

16.87 

15.74 

15.1 

14.59 

14.21 

18.4 


21.94 

20.54 

17.29 

16.09 

15.41 

14.9 

14.49 

18.9 


22.65 

21.26 

17.69 

16.41 

15.71 

15.18 

14.75 

19.4 


23.34 

21.98 

18.09 

16.71 

15.99 

15.44 

14.99 

19.9 



22.7 

18.49 

17.02 

16.27 

15.7 

15.22 

20.4 




18.92 

17.35 

16.55 

15.96 

15.45 

20.9 




19.36 

17.71 

16.85 

16.23 

15.69 

21.4 




19.82 

18.11 

17.17 

16.51 

15.93 

21.9 




20.29 

18.55 

17.5 

16.81 

16.18 

22.4 




20.77 

19.01 

17.84 

17.11 

16.43 

22.9 






18.19 

17.42 

16.69 

23.4 






18.53 

17.72 

16.95 

23.9 







18.03 

17.21 


Table 3. Empirical isochrones of the Pleiades in the J, H and Ks 

bands. More detailed version available on-line 


Yj 


J 

H 

Ks 

3.9 


4.09 


4.076 

4.4 

4.329 

4.517 

4.591 

4.532 

4.9 

4.831 

4.978 

5.054 

5.002 

5.4 

5.339 

5.444 

5.514 

5.47 

5.9 

5.862 

5.897 

5.957 

5.926 

6.4 

6.396 

6.329 

6.377 

6.35 

6.9 

6.946 

6.741 

6.768 

6.739 

7.4 

7.529 

7.119 

7.113 

7.086 

7.9 

8.118 

7.466 

7.416 

7.39 

8.4 

8.716 

7.775 

7.691 

7.65 

8.9 

9.315 

8.088 

7.962 

7.906 

9.4 

9.925 

8.423 

8.248 

8.185 

9.9 

10.56 

8.779 

8.561 

8.492 

10.4 


9.074 

8.871 

8.795 

10.9 


9.252 

9.143 

9.071 

11.4 



9.287 

9.239 


solar abundances as reference, see, e.g., [Allard et al.||2011) are 
all known to contribute to the differences between observations 
and model predictions. Such models should therefore be used 
with caution in the corresponding luminosity range. Whenever 
possible and until the models match the observations, empirical 
isochrones, such as the ones provided above should be preferred 
to their theoretical counterparts when looking for members. 


6. Present-day luminosity and mass functions of 
the Pleiades cluster 


photometry of the DANCe catalog was complemented by WISE 
( jCutri et al.|2012| l mid-infrared photometry, when available. 
Briefly, VOSA compares the observed SED to a grid of 


theoretical SEDs. The grid was made of BT-Settl models (AT 
lard et alI1|2012 1 in the range 1 500<Teff <4 000 K (in steps of 


100 K) and Kurucz models ( jCastelli et al.|[TW7] l in the range 
Teff >3 500 K. The overlap was defined on purpose to mitigate 
the effect of the discontinuity between the two sets of models 
over that range of effective temperatures. In that range, the best 
fit (in terms ofx^) was chosen between the two sets of models. 

VOSA offers the advantage of deriving Lboi using all the 
available photometric information rather than a subset of col¬ 
ors and luminosities. It estimates the bolometric luminosities by 
integrating the observed luminosities complemented by the best- 
fit model over the wavelength range not covered by the observa¬ 
tions. Eor most of the sample, at least 35% of the inferred bolo¬ 
metric luminosity corresponds to the observations. As such, it 
is less (but nevertheless still) sensitive to problems in a given 
photometric band pass, uncertain bolometric corrections, errors 
in the evolutionary models as described previously, variability, 
excesses, and some extent extinction. The Teff is computed us¬ 
ing the Lboi derived as described above and assuming a stellar 
radius. The latter depends mostly on the interior models, which 
have been relatively well calibrated for low mass stars. The Teff 
should therefore be only moderately affected by the deficiencies 
in the atmospheric models described in Section |fi!2] Using a sim¬ 
ilar method, pDahn et aLlp002] l estimated errors on the Teff on the 
order of 50 K, which is smaller than the sampling of the grid of 
models used in our study. 

The photometry was deredenned by assuming the canoni¬ 
cal mean extinction to the Pleiades of Ay =0.12 mag and a dis¬ 
tance of 136.2 pc (Melis et al. 2014|l was assumed to compute 


the bolometric luminosities. Photometric measurements below 
0.45 fim and beyond 12 pm are often subject to excesses related 
to activity and/or the presence of circumstellar debris disks and 
were systematically ignored for the fit. A visual inspection of 
the SEDs was performed to reject the remaining obviously prob¬ 
lematic photometric measurements. The quality and the quan¬ 
tity of the photometric measurements (~85% of the sources have 
more than 8 photometric measurements) resulted in a robust fit 
for most sources. Only six sources had fewer than the minimum 
three photometric measurements required to perform a meaning¬ 
ful fit, and are discarded for the rest of the analysis. The DANCe 
Pleiades Release 2 online table includes the estimated Teff and 
Lboi for the 2 010 sources selected as members in this study. 


6.2. Completeness and contamination of the sample 

Prior to deriving a luminosity function, we assess the complete¬ 
ness of our analysis. We identify several sources of incomplete¬ 
ness and contamination. As described in jSarro et al.| ( |2014l l, the 
75% selection threshold on the membership probability leads to 
an estimated contamination rate of ~7% at most. The number of 
missed members is estimated to reach 4% at most. These num¬ 
bers do not take real Pleiades members with properties differing 
significantly from their siblings into account; 


6.1. Luminosities and effective temperatures of the Pleiades 
members 

We estimated the effective temperatures (hereafter Tetr) and 
bolometric luminosity (Lboi) of all Pleiades candidate members 
using the virtual observatory SED analyzer (hereafter VOSA, 
[Bayo et al.|2008|l. Eor that purpose, the optical and near-infrared 


- High order (N >3) unresolved multiple systems can appear 
significantly brighter than the single star and equal-mass bi¬ 
nary sequences used for the calculations and hence be given 
lower probabilities. Eigure shows that several candidate 
members passed our selection criterion in spite of being over- 
luminous and could be high order multiple systems. 
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Fig. 7. Comparison of the observations of the Pleiades sequence (black dots) and models from | Allard] ( |2014) ) (red) and ||Bressan et al.| ( |2012| l (light 
green) for 120 Myr in various color magnitude diagrams. 
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- The orbital motion of long period (decades or centuries) un¬ 
resolved binaries can affect their apparent proper motion sig¬ 
nificantly and result in lower membership probabilities. 

- White dwarfs have colors and luminosities inconsistent with 
the empirical stellar main sequences used to compute the 
membership probabilities and were missed. A future study 
will look for them in the DANCe catalog. 


Additionally, because the DANCe catalog is based on 
archival observations of very diverse origins, the sensitivity - 
hence completeness - is not homogeneous over the entire area 
covered by the survey. Based on the spatial coverage and depth 
of the datasets used to derive the DANCe catalog, we identify 
a region of the survey where the completeness is believed to 
be homogeneous. This area is defined by our own recent CFHT 
MegaCam f-band observations and covers a roughly square re¬ 
gion of ~3x3° centered on the Seven Sisters (see Fig. 1 of |Bouy 


|et al.||20T3] l. To perform a meaningful study of the luminosity 
and mass functions, we therefore limit the following analysis to 
a circular region around the cluster’s center (at RA=03:46:48 and 
Dec=24:10:17 J2000) and with a radius of 3°. A total of 1 378 
members fall within this region. 

[Pinfield et al.| ( |1998| ) reported a significant mass segregation 
that might affect the present study of the mass and luminosity 
function. We nevertheless note that the 3° radius corresponds 
to more than double the size of the core radius for low mass 
stars and almost seven times the core radius for high mass stars 
estimated by Pinfield et al. ( 1998|l. 


6.3. System’s present-day luminosity function 

Figure shows the bolometric and absolute Ks-band system’s 
luminosity function assuming a distance of 136.2 pc. We chose 
to represent the distribution of absolute Ks-band magnitudes be¬ 
cause all but 8 of the 1 378 members falling inside the region 
defined above have a Ks-band measurement.The DANCe cata¬ 
log is believed to be sensitive up to Ks~20 mag but complete 
up to Ks~17 mag (or M^s =11.3 mag). A kernel density esti¬ 
mate (KDE) was preferred over histograms to build a continu¬ 
ous present-day luminosity function (PDLF). The KDE was ob¬ 
tained using a bandwidth of 0.3 mag for the DANCe sample and 
0.8 mag for the Tycho-2 sample and an Epanechnikov kernel. 
Uncertainties were estimated using bootstrapping and comput¬ 
ing the 95% confidence intervals (defined between quantile 2.5% 
and 97.5%) over 1 000 replicates. These uncertainties do not in¬ 
clude the contamination and missed members mentioned above 
or the photometric uncertainties. 

An interesting feature is observed beyond Mks >9.8 mag 
(Ks>15.5 mag) and Mboi >12.5 mag, where the luminosity func¬ 
tion becomes almost flat. These luminosities correspond roughly 
to 0.04~0.05 Mq according to the models of [Allard (2014|l as¬ 


suming an age of 120 Myr ( [Stauffer et al. 1998| Barrado y Navas 
cues et aL]|2015 i. This significant change in the slope suggests 


that different or additional formation mechanisms must be at 
work for high and low mass brown dwarfs. 

We note that the small drop at the edges of the complete do¬ 
main is interpreted as a boundary effect of the kernel density es¬ 
timate and is not physical. The luminosity function is available 
in the online Table (Tj 


6.4. System’s present-day mass function 

Deriving a mass function requires transforming luminosities to 
masses. These transformations are usually based on poorly or 


uncalibrated theoretical mass-luminosity relationships. In Sec¬ 
tion we found that the models poorly match the observed 
colors and luminosities, especially at the very low mass end. 
The mass-luminosity relationship in the corresponding luminos¬ 
ity range must therefore be considered with caution, and any 
mass estimate should be regarded as strongly model dependent 
and tentative. With this important limitation in mind, we derive 
the Pleiades present-day mass function for system (PDME) and 
compare it to recent predictions and models of the system’s ini¬ 
tial mass function (IMP). 

The values of Lboi estimated with VOSA were transformed 
into masses using the mass-luminosity relationship given by the 
120 Myr BT-Settl models for objects with Tgff <4000 K and 
the 120 Myr Bressan et al.| ( |201^ models for objects hotter than 
Teff >4000 K. According to these calculations, the faintest can¬ 
didate members have masses of only ~20 Mjup, and the most 
massive member is Alcyone with 9.4 Mq. The estimated masses 
are also included in the online DANCe Pleiades Release 2 table. 

Pigure[9]shows the distribution of masses obtained. The KDE 
was obtained using a bandwidth of 0.1 (in logio(m)) for the 
DANCe and Tycho-2 samples and an Epanechnikov kernel in 
both cases. As for the PDLP, uncertainties were estimated us¬ 
ing bootstrapping and computing the 95% confidence intervals 
(defined between quantiles 2.5% and 97.5%) over 1 000 repli¬ 
cates, and they do not include the contamination and missed 
members mentioned in Sectio n|6.2| Table|2]gives the results, and 
Pig.|9 compares to Chabrier ( |2005| l model and jThies & Krou^ 
P0()7 I prediction normalized to the same total mass over the 
range 0.03-0.6 Mq. Here as well the small drop at the edges of 
the complete domain is interpreted as a boundary effect of the 
kernel density estimate so is not physical. 


6.5. Comparison to theoretical predictions 


Unlike what was reported by previous authors (see, e.g., Moraux 
et al.|2003) Kroupa et al.|2013|l, the slope abov e 1 Mq matches 
the a=2.3 value predicted by Thies & Kroupa P007| l relatively 
well. But as mentioned above, the study of membership by 


Sarro et al. (2014 1 used for the present analysis did not include 
white dwarfs. At least two white dwarf candidate members of 
the Pleiades were reported in the literature (Dobbie et al.|2006 


Luyten & Herbig|1960| l, one of them ultra-massive. If confirmed, 
these two members could significantly change the slope of the 
PDME at the high mass end. Additionally, early dynamical inter¬ 
actions in multiple systems probably ejected a substantial num¬ 
ber of massive stars. The missing white dwarfs and ejected mem¬ 


bers therefore affect the comparison with the model of Chabrier 
P005)l and Thies & Kroupaj (|2007|l, and suggest that in spite 


of the apparently better match of Chabrierj ( [2005 | l model with 
our observations, [Thies & Krou^ ( [2007[ l prediction is probably 
closer to the real distribution (in that mass range). 

Overall, jChabrierl ( [2005[ l model and [Thies & Krou^ ( [2007[ l 
match the PDME in the range 0.02<M<0.6 Mq relatively well. 
Both predict too many very low mass stars and brown dwarfs, 
a discrepancy that has been repeatedly reported in the Pleiades 
and other clusters and could be related to uncertainties in the 
mass-luminosity relationship used to build the observed PDME 
Interestingly, the PDME becomes flat at the very end, between 
20 Mjup and 50 Mjup, a feature already suggested by other au¬ 
thors and now confirmed by our study and the consequence of 
the flattening observed in the luminosity function mentioned 
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Fig. 8. Top panel: Observed present-day bolometric luminosity function for systems within the central 3° of the Pleiades cluster. Bottom panel: 
Observed present-day absolute Ks-band luminosity function for systems within the central 3° of the Pleiades cluster. Domains where the survey 
is incomplete are indicated by shaded areas. 


above. One possible explanation was that several formation 
mechanisms lead to the formation of brown dwarfs. While the 
majority of brown dwarfs probably form like their more mas¬ 
sive stellar counterparts from the fragmentation and collapse of 
a molecular cloud, a significant number might be formed as com¬ 
panions in the disk of a more massive system and subsequently 
ejected by dynamical interactions. 


As illustrated in Fig. the observations suggest dividing 
the PDMF into four domains. Each domain can be fitted with a 
power law with the following indices; 

1. logio(m) >0.2 : ai- 2.56 (missing the white dwarfs); 

2. -0.7< logio(m) <-0.22 : a 2 = 0.23; 

3. -1.3< logiolm) <-0.7 : q' 3 = -1.22; 

4. -1.45</ogfio(m) <-1.3 : a4= 0.13. 
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Fig. 9. Observed present-day system mass function of the Pleiades within the central 3° radius (red). |Chabrier| |2005| l model for the galactic disk 
mass function (normalized to the same total mass over the range 0.03-0.6 Mq) and |Thies & Krou^j2007^ prediction for the Pleiades cluster are 
overplotted in blue and green, respectively, and normalized to the same total mass over the range 0.03-0.6 Mq. The “gap” between the Tycho-2 
and DANCe analyses and the estimated limit of completeness of the DANCe survey are indicated by shaded areas. 


6.6. The mass function in time: comparison to other dusters 


Figure o compares the results of our study to the PDMF re¬ 
ported in the young ~1 Myr old Trapezium cluster and in the 
625 Myr old Hyades cluster. The Hyades PDMF was computed 
as follows. The stellar PDMF was computed using the sample of 
members and masses reported in |Goldman et al.| ( |2013| l and the 
substellar PDMF using the sample of brown dwarfs and masses 
reported in Bouvier et al. ( 2008| l. In both cases a KDE with an 
Epanechnikov kernel with a bandwidth of 0.1 (in log\o{m)) was 
used. The two studies covered very different areas. The substel¬ 
lar PDMF obtained from the Bouvier et al. (2008i sample was 
therefore scaled to match the integrated mass of the |Goldman| 
|et al.| ( |2()13[ ) sample over the mass range 0.1-0.15 Mq where 
both studies are expected to be complete and overlap. 

The PDMF of the Trapezium cluster was estimated by trans¬ 


forming the /T-band luminosity function reported by Muench 
et al.| (|2002 1 into a mass function using the BT-Settl models at 


1 Myr up to 1 Mq, and |Bressan et al.| ( |2012[ ) at 1 Myr models 
above 1 Mq. 


The Trapezium and Hyades clusters were not chosen ran¬ 
domly for the comparison. They are similar to the Pleiades in 
many aspects. Their total masses and numbers of members are 
similar, and various studies have derived their mass functions 
over a similar mass domain, making the comparison meaningful 
and interesting. Additionally, the Pleiades has been described as 
a snapshot of the Trapezium’s future (|Kroupa et al.||20()T]l and 


of the Hyades’ past (Bouvier et al. 2008) 1. Figure [TT| therefore 
illustrates the evolution of the mass function in time. 


In this context and as already known (see, e.g., IMoraux &| 
|Bouvier|2012| and references therein). Figure [TT|shows that the 
proportion of very low mass stars and brown dwarfs is clearly 
decreasing over time and that the peak of the mass function is 
shifting toward higher masses from 0.2 Mg at the age of the 
Trapezium and of the Pleiades to 0.5 ~1 Mg by the age of the 
Hyades. 

The flattening of the Pleiades PDMF below 0.05~0.06 Mg is 
observed in the Trapezium PDMF as well, reinforcing the idea 
that the overabundance of low mass brown dwarfs is due to ad¬ 
ditional formation mechanisms - including photo-evaporation, 
disk instabilities, and ejections - rather than late (>1 Myr) dy¬ 
namical evolution. A similar flattening might be present in the 
Hyades PDMF around 0.075 Mg, but the very large uncertain¬ 
ties due to the very-small-number statistics in that mass range 
prevent any Arm conclusion and meaningful comparison. 

In any case, and within the large uncertainties of each PDMF 
at the very low mass end, the level of overabundance around 
0.05 Mg seems to be decreasing with time, since it is rather large 
at the age of the Trapezium and the Pleiades and possibly shifted 
toward higher masses (if present) at the age of the Hyades. This 
would be consistent with the classical picture in which a signifi¬ 
cant number of very low mass stars and brown dwarfs are ejected 
before the first million years, and subsequently evaporate as the 
cluster becomes more and more relaxed. 
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Fig. 10. Observed present-day system mass function of the Pleiades within the central 3° radius (red). Best fit power laws in the form m “ are 
overplotted in blue, and their index is indicated. The “gap” between the Tycho-2 and DANCe analyses is indicated by shaded areas. 


At the high mass end, the comparison is not straightforward. 
Both the Hyades and Pleiades samples miss the white dwarfs. 
Including the white dwarfs would certainly change the compar¬ 
ison. The Trapezium includes several massive stars. By the age 
of the Pleiades, the most massive of them (which fall beyond 
Fig. [TT|coverage) have since long exploded as supernovae . Some 
massive stars might also have been ejected from their parent mul¬ 
tiple stellar system. In any case, the lack of statistics at the high 
mass end means that the comparison of the mass function over 
the corresponding domain is often not statistically significant. 


7. Conclusions 


We have presented a study of the Pleiades cluster’s properties. 
We combined the results of the DANCe survey (Satro et al. 2014[ 
updated with new photometry and improved membership prob¬ 
abilities) with a reanalysis of the Tycho-2 catalog over the same 
area. Our analysis allows us to: 


1. identify 2 010 candidate members in the DANCe survey; 

2. identify 207 candidate members in the Tycho-2 catalog, in¬ 
cluding 87 new members; 

3. provide the largest and most complete census of the cluster. 
After removing matched entries in the DANCe and Tycho- 
2 samples (110 sources) and adding Electra and Alcyone 
(missed in both surveys), the final sample of Pleiades mem¬ 
bers adds to 2 109 members, including 812 previously un¬ 
known members; 

4. estimate contamination levels between 35% and 61% in 
previous surveys from the literature searching for Pleiades 
brown dwarfs; 
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5. provide accurate empirical isochrones in six broad band op¬ 
tical and near-infrared filters ir,i,Y,J,H,Ks)\ 

6. show that theoretical isochrones poorly match the empirical 
sequence in most color-magnitude diagrams made of optical 
and near-infrared colors and luminosities; 

7. observe that the luminosity function becomes flat beyond 
Ks ^15.5 mag, corresponding to ~0.04 - 0.05 Mq and sug¬ 
gesting that different or additional formation mechanisms 
must be at work at the very low mass end; 

8. show that the overall shape of the theoretical models or pre¬ 
dictions of the system initial mass function is in reason¬ 
able agreement with the observed system’s present-day mass 
function (assuming that the theoretical mass-luminosity re¬ 
lationship is correct). They predict too many brown dwarfs, 
but the uncertainty on the mass-luminosity relationships pre¬ 
vents us from drawing any Arm conclusion about this differ¬ 
ence. They do not predict the flattening down to 0.05 Mq 
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Fig. 11. Observed present-day system mass function of the Pleiades within the central 3° radius (red), compared to the measurements reported for 
the Trapezium cluster by|Muench et al.|p002| t and |Hillenbrand & Carpent^ ( |20Q0| > (blue) and the Hyades (orange) from |Goldman et ^j2Q13f 
(down to 0.15 Mo) and |Bouvier et al■H2008^ in the range 0.05<M<0.15 MO. 
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Table 1. Comparison with previous Pleiades brown dwarf surveys 



Moraux+2003 

Biham-(-2UU6 

Casewell-l-2007 

ZO+2024 

RA 

Dec 

San'o-i-2014 

DANCe 

proba 

Comments 

CFHT-PLIZ-1 




57.774916 

24.604747 

5277684 

J03510596+2436170 

0.1 


CFHT-PLIZ-2 


PLZJ-78 


58.846125 

24.818058 

5296873 

103552307-1-2449051 

1.0 


CFHT-PLIZ-3 

BRBll 



58.028 

24.266877 

5172848 

J03520671-P2416006 

1.0 


CFHT-PLIZ-4 




55.4205 

25.906388 

5384327 

J03414090-P2554242 

1.0 


CFHT-PLIZ-5 




58.408166 

26.038797 




(1) 

CFHT-PLIZ-6 

BRB9 

PLZJ-9 


58.479583 

23.393447 

5061235 

J03535512-f2323361 

1.0 


CFHT-PLIZ-7 




57.050541 

25.907888 

5393326 

J03481212-t-2554283 

0.046 


CFHT-PLIZ-8 




55.75075 

24.731147 

5224304 

J03430017-P2443525 

1.0 


CFHT-PLIZ-9 




56.146625 

25.228427 

5236418 

103443516-1-2513429 

1.0 


CFHT-PLIZ-10 




57.937375 

23.444297 

5050329 

103514494-f2326395 

1.0 


CFHT-PLIZ-12 

BRB14 



57.856708 

23.755877 

5163704 

103512557-P2345214 

1.0 


CFHT-PLIZ-13 




58.768333 

26.263699 

4292417 

103550436-1-2615503 

1.0 


CFHT-PLIZ-14 




58.384958 

26.116999 




(1) 

CFHT-PLIZ-15 




58.077666 

24.074558 

5170275 

103521864-P2404280 

1.0 


CFHT-PLIZ-16 




55.917875 

24.50315 

5227533 

103434030-P2430114 

1.0 


CFHT-PLIZ-17 




57.861208 

23.502958 

5051897 

103512665-^2330108 

0.996 


CFHT-PLIZ-18 




58.503999 

24.914697 




(1) 

CFHT-PLIZ-19 




59.068208 

23.914288 

5181049 

103561638-1-2354514 

0.999 


CFHT-PLIZ-20 

BRB15 

PLZJ-11 


58.522375 

23.566519 

5169202 

103540534-1-2333595 

1.0 


CFHT-PLIZ-21 




58.865249 

25.827977 

5406704 

103552762-f2549407 

0.977 


CFHT-PLIZ-22 




57.969624 

26.875599 

5535893 

103515262-P2652319 

0.0 


CFHT-PLIZ-23 




57.8895 

24.170599 

5161880 

103513347-t-2410142 

0.0 


CFHT-PLIZ-24 




56.848666 

26.016597 




(1) 

CFHT-PLIZ-26 




56.202749 

25.654866 

5390235 

103444868-P2539175 

0.0 


CFHT-PLIZ-28 

BRB18 



58.558458 

23.297608 

5054020 

103541406-1-2317521 

1.0 


CFHT-PLIZ-29 




57.438708 

26.847188 

5518534 

103494524-P2650503 

0.0 


CFHT-PLIZ-30 




57.941666 

26.827058 

4177219 

103514595-1-2649382 

0.0 


CFHT-PLIZ-31 


PLZJ-2I 


57.948541 

24.66653 

5285380 

103514765-^2439591 

1.0 


CFHT-PLIZ-32 




57.564458 

26.580908 

5517035 

103501543-^2634515 

0.0 


CFHT-PLIZ-33 




57.686166 

26.702599 




(1) 

CFHT-PLIZ-34 




58.510666 

24.673908 

5295295 

103540255-^2440259 

0.0 


CFHT-PLIZ-35 

BRB20 



58.163208 

24.775008 

5290268 

103523915-^2446296 

1.0 


CFHT-PLIZ-36 




58.65975 

23.633508 

5169338 

103543837-f2338011 

0.0 


CFHT-PLIZ-37 




58.914875 

24.214477 




(1) 

CFHT-PLIZ-3 8 




56.477875 

26.504047 

5512940 

103455470-1-2630149 

0.0 


CFHT-PLIZ-39 




58.417916 

26.271708 

4181633 

103534029-1-2616198 

0.0 


CFHT-PLIZ-40 




57.455416 

26.565608 

5518403 

103494925-f2633563 

0.0 



BRBl 



57.964999 

23.580611 

5161183 

103515155-1-2334492 

1.0 



BRB2 



58.184583 

23.904222 

5171639 

103524428-1-2354151 

0.999 



BRB3 



58.032916 

23.987388 




(1) 


BRB4 

PLZJ-29 


56.096666 

25.645749 

5385058 

103442323-^2538451 

0.997 



BRB5 



58.215833 

23.563583 

5174064 

103525179-f2333482 

0.999 



BRB6 



58.29 

23.563416 

5173716 

103530963-1-2333479 

0.999 



BRB7 



57.18625 

24.622972 

5279607 

103484469-^2437236 

1.0 



BRB8 



58.024166 

24.292138 

5172897 

103520582-f2417312 

1.0 



BRBIO 



57.312916 

24.606222 

5277273 

103491512-P2436225 

1.0 



BRB12 



56.9125 

24.606138 

5238416 

103473901-P2436223 

1.0 



BRBl 3 



58.802083 

23.293888 

5058757 

103551261-f2317373 

0.0 



BRB16 



57.130833 

24.577138 

5277223 

103483153-^2434372 

1.0 



BRB17 



58.534624 

23.909277 




(1) 


BRB19 



58.714541 

23.753388 

5177022 

103545147-P2345122 

0.0 



BRB21 

PLZJ-4 


58.542791 

23.694527 

5169447 

103541027-P2341402 

1.0 



BRB22 

PLZJ-61 


56.130291 

25.587527 

5390123 

103443127-f2535149 

1.0 



BRB23 



57.664708 

25.048472 

5278239 

103503955-^2502545 

1.0 



BRB24 



57.081875 

24.869305 




(1) 


BRB25 



58.248875 

23.865555 

5176839 

103525971-P2351557 

0.0 



BRB26 



57.234041 

25.161972 




(1) 


BRB27 

PLZJ-32 


56.113499 

25.744972 

5388791 

103442728-1-2544419 

0.996 



BRB28 

PLZJ-37 


58.228833 

24.621833 

5290218 

103525491-^2437181 

0.999 



BRB29 



58.505958 

23.832805 

5171913 

103540143-f2349575 

0.998 



BRB30 



56.987333 

24.482888 




(1) 


BRB31 



57.858041 

24.258944 




(1) 


BRB32 



57.253583 

24.694027 




(1) 


BRB33 



57.833958 

23.755111 

2535230 

103512013-1-2345178 

0.0 



BRB34 



57.193958 

24.750888 




(1) 



PLZJ-56 


56.222958 

25.605405 

5390145 

103445352-1-2536194 

0.0 




PLZJ-45 


58.2425 

24.292102 







PLZJ-50 


55.983249 

25.607069 

4055668 

103435600-1-2536253 

1.0 




PLZJ-60 


56.134666 

25.421683 

5387838 

103443232-1-2525181 

1.0 




PLZJ-46 


58.238333 

24.266944 







PLZJ-77 


58.541833 

23.297855 







PLZJ-10 


58.129958 

24.774891 







PLZJ-23 


57.972416 

24.636697 







PLZJ-93 


58.804166 

24.604388 

2629741 

103551301-1-2436158 

0.0 




PLZJ-323 


55.980291 

25.723944 

2646887 

103435527-1-2543262 

0.0 




PLZJ-721 


58.779749 

24.956205 

2618690 

103550717-1-2457230 

0.044 




PLZJ-235 


58.135708 

24.743511 







PLZJ-112 


58.330708 

24.89218 

2836575 

103531938-1-2453317 

0.0 




PLZJ-100 


56.829958 

25.848138 








PNMOl 

57.174958 

22.916166 

5046135 

103484198-1-2254583 

0.0 





PNM02 

55.461541 

25.276361 

5231515 

1034150764-2516351 

0.0 





PNM03 

57.02925 

22.944361 

5006623 

1034807024-2256396 

0.0 





PNM04 

57.130666 

25.715861 

5399403 

1034831354-2542571 

0.0 





PNM05 

57.013291 

23.201611 

5005906 

1034803194-2312057 

0.0 
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PNM06 

55.525875 

25.075277 

5231041 

J03420620+2504310 

0.0 


PNM07 

55.466583 

25.216805 

5231375 

J03415198+2513005 

0.0 


PNM08 

57.134416 

25.715833 

5399402 

J03483226+2542569 

0.0 


PNM09 

55.781 

25.266638 

5225601 

J03430744+2515598 

0.0 


PNMIO 

56.991666 

25.820138 

5393154 

J03475799+2549125 

0.0 


PNMll 

57.016458 

22.876999 

5007559 

J03480394+2252372 

0.0 


PNM12 

55.69625 

25.199194 

5225443 

J03424710+2511573 

0.0 


PNM13 

55.789624 

25.268361 

4006442 

J03430951+2516060 

0.0 


PNM14 

57.001333 

25.690083 

4025181 

J03480033+2541242 

0.0 


PNM15 

57.162833 

22.887749 

4000699 

J03483909+2253158 

0.0 


Calar Pleiades 8 

55.892458 

25.037944 

5226398 

J03433419+2502166 

0.96 


Calar Pleiades 9 

55.512875 

24.061416 

5144919 

J03420308+2403413 

0.0 


Calar Pleiades 10 

56.558583 

23.365666 

5008540 

J03461406+2321565 

1.0 


Calar Pleiades 11 

55.575083 

24.919444 

5228413 

J03421802+2455099 

1.0 


Calar Pleiades 12 

55.70075 

24.066999 

5139930 

J03424817+2404012 

1.0 


Calar Pleiades 13 

55.433083 

25.058472 

5228741 

J03414395+2503305 

1.0 


Calar Pleiades 14 

56.460958 

23.743555 

5155777 

J03455062+2344369 

1.0 


Calar Pleiades 15 

55.627458 

25.04425 

5228710 

J03423058+2502392 

1.0 


Calar Pleiades 16 

56.944916 

25.587944 

5394991 

J03474678+2535166 

1.0 


Calar Pleiades 17 

56.298874 

23.695444 

5148223 

J03451172+2341436 

1.0 


Calar Pleiades 18 

56.388749 

23.576194 

5153508 

J03453330+2334342 

1.0 


Calar Pleiades 19 

57.065208 

25.835777 

5393192 

J03481565+2550089 

0.999 


Calar Pleiades 20 

56.493666 

23.698277 

5153687 

J03455848+2341539 

0.759 


Calar Pleiades 21 

56.7515 

23.365777 





Calar Pleiades 22 

56.501625 

23.326083 





Calar Pleiades 23 

55.729708 

24.150722 




(2) 

Calar Pleiades 24 

56.702791 

23.480194 





Calar Pleiades 25 

57.110916 

22.864972 





Calar Pleiades 26 

56.317041 

23.585361 

5148179 

J03451608+2335088 

0.0 

misidentification 

Calar Pleiades 27 

56.550541 

23.304361 

3387884 

J03461212+2318156 

0.0 


Calar Pleiades 28 

56.657875 

23.424249 




(2) 

Calar Pleiades 29 

55.72575 

24.136527 




(2) 

Calar Pleiades 30 

57.006208 

25.613666 





Calar Pleiades 31 

56.427041 

23.730388 




(3) 

R1 

55.753 

25.084611 

3446081 

J03430072+2505046 

0.0 


R2 

57.005416 

23.136972 




(2) 

R3 

57.117416 

23.209166 




(2) 

R4 

57.456708 

22.898749 

3430893 

J03494961+2253556 

0.0020 


R5 

57.076333 

25.456305 

2770428 

J03481832+2527227 

0.0 



Notes. (1) No detection even though the survey was sensitive enough to detect the source; (2) i-band luminosity inconsistent with membership; (3) 
r-band luminosity inconsistent with membership; 
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Appendix A: DANCe-Pleiades Catalogue Data 
Release 2 

In this new release, the DANCe-Pleiades catalogue presented in 
|Bouy et ah] ( |2013| l and ( |Sarro et al.|2014[ ) has been improved in 
the following ways: 

- addition of APASS gri photometry when no measurement 
was available; 

- addition of F-band photometry for 40 bright candidate mem¬ 
bers from the list of |Stauffer et al.| ( [2007| l; 

- recalculation of the membership probability using the new 


photometry following the method described in (Sarro et al. 

[20T4l l. 

Complementing the previous release at the bright end using 
the APASS and our own F-band observations was particularly 
useful for the membership analysis. The selection process de¬ 
scribed in [Sarro et al.| ( |2014] l made use of the ri and F JHK 
photometry. The JHK photometry comes from a combination 
of 2MASS and UKIDSS and covers the entire dynamic range up 
to ~19 mag. On the other hand, the riY photometry was limited 
to ~ 13 mag, and was preventing us from defining a training set of 
members with complete photometric coverage for the brightest 
sources. 

To circumvent this problem, we gathered gri photometry 
from the APASS survey. The APASS DR7 release covers ~91% 
of the sky and 100% of the area encompassed by the DANCe- 
Pleiades survey in the Johnson B and V and Sloan gri biters. 
Since no Johnson photometry was available in the DANCe- 
Pleiades catalog, it was ignored and we only included the gri 
photometry. The APASS survey is valid from about the 10'^ to 
17'* magnitudes. Brighter sources are not as reliable, but we nev¬ 
ertheless chose to report the corresponding photometry in our 
catalog, keeping this limitation in mind. 

No shallow F-band survey was available in the Pleiades at 
the time of our analysis. Service time was granted to us at the 
William Herschel Telescope (Proposal SW2013b27, PI. Bouy) 
with its LIRIS near-infrared camera ( jAcosta-Pulido et al.|2002| l. 
Shallow F-band images of a sample of 40 bright candidate mem¬ 
bers of [Stauffer et al.| ( [2007| l with available riJHK photometry 
and randomly selected to cover the luminosity range between 
8-14 mag were obtained on 23 December 2013 and 3 January 
2014. Sets of bve individual exposures of 0.9 s were obtained 
at bve positions around the target. The ambient condition was 
good, and the photometric zeropoints were derived by compari¬ 
son of the faint sources present in the images with their UKIDSS 
counterparts. The individual raw images were processed using 
an updated version of Alambic (Vandame 2002[ ), a software 
suite developed and optimized for the processing of optical and 
near-infrared imagers, and conbgured for LIRIS. Alambic in¬ 
cludes standard processing procedures, such as dark subtraction 
for each individual readout ports, bat-beld correction, bad pixel 
masking, destripping, and background subtraction. 

The analysis presented in Sarro et al. ( [2014 1 was then re¬ 
peated including the new photometry. It allowed us to extend the 
training set, hence the bnal calculation of membership probabil¬ 
ities, at the high luminosity end up to ~9 mag. As a result, a total 
of 97 members originally missed by our analysis were recovered 
at the bright end. Table gives the astrometry, photometry and 
membership probability for the second release of the DANCe 
catalogue presented in this study. 

Figure [ATT] shows the disttibution of apparent magnitudes in 
the i and Ks bands for the entire sample. These histograms pro¬ 
vide a useful estimate of the completeness and sensitivity limits 



Fig. A.l. Distribution of apparent magnitudes for all the sources in the 
survey in the /-band (blue) and Ki-band (red). 


of the survey. In the /-band, the survey is sensitive up to 25 mag 
and complete up to ~23 mag. In the Ki-band, the survey is sen¬ 
sitive up to ~20 mag and complete up to ~18 mag. A small ex¬ 
cess appears around i ~15 mag and Ks ~14 mag. In the /-band, 
it corresponds to the overlap between APASS and the DANCe 
photometry, with the APASS survey encompassing a larger area. 
In the Ks-band, it corresponds to the overlap between 2MASS 
and UKIDSS, which have different spatial resolutions. 


Appendix B: Selection of Pleiades members In the 
Tycho-2 catalog 

Although the Tycho-2 catalog has been extensively used to 
search for Pleiades members, we decided to perform our analysis 
using the method developed by [Sarro et al.[ ( [20T4| ). We retrieved 
the Tycho-2 catalog within the same area as the DANCe sur¬ 
vey and complemented it with APASS, 2MASS, CMC-14, and 
WISE photometry. The best set of variables (in a statistical sense, 
see [Sarro et al.|2014[ for more details) was found to include the 
proper motions and the Bj, Vt, r, J, H, and K magnitudes. The 
training set was made of 175 candidate members from the list of 
[Stauffer et al.[ ( 2007) with a counterpart in the Tycho-2 catalog. 
Simulations of the completeness and contamination (see [Sarro[ 
[et al.|201^ for more details) suggest to use a threshold of 48% 
of probability to select members in the hnal catalog. Figure jBTTj 
shows a Bj vs Bt - Vr and vector point diagrams for the 5 537 
sources included in this study. With this criterion, a total of 207 
candidate members are selected. We find that: 


83 of the 207 candidate members were not present in [Stauffer 
[et al.[([2007ll list, and 


19 cahdidate members out of the original 175 of Stauffer 
[et al.[ ( |2007) l list with a Tycho-2 counterpart are not selected 
in our analysis. Figure jBTTj shows that a significant fraction 
of them are rejected based on inconsistent proper motions, 
and a few of them are only based on inconsistent colors and 
luminosities. 

Tablej^gives the astrometry, photometry and membership prob¬ 
ability for the 5 537 Tycho-2 sources included in this study. 
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Fig. B.l. Upper panel: Bj vs Bj - Vr color-magnitude diagram of the 
Tycho-2 catalog in the area of the DANCe survey (black dots). New 
candidate members (with membership probability greater than of equal 
to 48%) are represented by red dots. Candidate members from |Stauffer| 
|et al.| j2007^ that do not pass the selection criterion are represented by 
green squares. Lower panel: vector point diagram of the Tycho-2 cata¬ 
log in the area of the DANCe survey. The color and symbol codes are 
the same. 


Article number, page 18 of 


20 













A(&A-AA-2014-25019, Online Material p 19 


Table 1. Present day system bolometric luminosity function of the Pleiades (normalized) 


Mboi (MO) 


Density 


0.0018+^)!^,-'* 
0 . 0022+0 
0.0025+ 
0.0036+ 
0.0066+° 
0.011+°-'^ 
0.0167+ 
0.0215+ 
0.027+^ 
0.0277+'®-°'^ 
0.0831+° 

0 . 1088 + 0^27 
0 - 1796 !± 
0 . 2668 !°gi 

0 9447+0 0634 

2 ri'tnc 

0.1234+; 
0.0615+; 
0.0388^; 
0.0397+: 


-3.14 

-2.26 

-1.38 

-0.5 

0.38 

1.26 

2.14 

3.02 

3.89 

4.77 

7.21 

8.04 

8.86 

9.69 

10.51 

11.33 

12.16 

12.98 

13.8 


0137 


Table 2. Present day system mass function of the Pleiades (normalized) assuming an age of 120 Myr, a distance of 136.2 pc and t 
extinction of Ay =0.12 mag 


logmjm) ^(m) 


-1.55 

0.18+°° 

-1.41 

0.24+°- 

-1.26 

0.25+1 

-1.11 

0.34!°- 

-0.96 

0.54+1 

-0.82 

0.85+|§ 

-0.67 

1.08^1 

-0.52 

1.05^1 

-0.38 

0.96+1 

-0.23 

0.81+1 

0.13 

0.29+°+ 

0.19 

0.24+°- 

0.25 

0.19+1 

0.31 

0.14+|[ 

0.37 

0.1+°| 

0.43 

0.07+°+ 

0.5 

0.05+1 

0.56 

0.03+°+ 

0.62 

0.02!°+ 

0.68 

0 


mean 
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Table 3. DANCe catalogue of the Pleiades (second release) 


DANCe 

Sarro+2014 

NPOS_TOTAL 

NPOS_ 

_OK 

RA 

Dec 

err_RA 


J035008.43+253255.5 

5403154 

44 


44 

57.5351086587683 

25.5487606105672 

3.3705095E-7 

3.3' 

J034623.03+243617.8 

5237969 

80 


80 

56.5959755649232 

24.6049563409070 

2.4088035E-7 

2.4( 

J034524.79+242045.1 

5157061 

75 


75 

56.3532988755393 

24.3458727345082 

2.5063028E-7 

2.: 

J035254.91+243718.1 

5290218 

46 


46 

58.2288026911685 

24.6216882239095 

1.2403241E-6 

1.2‘ 

J034133.72+250648.5 

5224174 

41 


41 

55.3904898811303 

25.1134754192527 

2.0917902E-7 

2.( 

J034043.20+224953.8 

4991637 

45 


45 

55.1799876390125 

22.8316197776090 

1.6703784E-7 

1.6' 

J034809.22+235840.5 

5167474 

146 


146 

57.0384110265741 

23.9779110883695 

2.0908E-7 

2.0' 

J034744.44+243655.6 

5238429 

108 


108 

56.9351582236966 

24.6154514381828 

5.5228713E-7 

5.5: 

J034505.00+234606.4 

5150226 

83 


83 

56.2708149165531 

23.7684469977015 

2.404183E-7 

2.4( 

J034651.44+240616.0 

5149477 

94 


94 

56.7143228114129 

24.1044497268851 

2.8111455E-7 

2.8 

J034336.58+231233.9 

4997214 

49 


49 

55.9024017319216 

23.2094137491369 

3.154203E-7 

3. 

J034342.89+255137.0 

5383859 

46 


46 

55.9287136644351 

25.8602804825768 

2.2304917E-7 

2.2: 

J034411.28+245234.1 

5227884 

61 


61 

56.0470181374588 

24.8761386932101 

4.078212E-7 

4.( 

J035054.65+242155.9 

5157473 

94 


94 

57.7276985353782 

24.3655351189338 

1.2143116E-7 

1.2 

J035244.28+235415.1 

5171639 

42 


42 

58.1845187631600 

23.9042076976708 

1.9325248E-7 

1.9: 

J034227.31+223424.5 

4935291 

39 


39 

55.6138046074025 

22.5734844985916 

5.0468094E-7 

5.0‘ 

J034427.50+241416.9 

5152272 

73 


73 

56.1145974564091 

24.2380192856017 

2.511189E-7 

2.5 

J034024.19+243504.0 

4134727 

56 


56 

55.1008096693649 

24.5844530665781 

1.6823185E-7 

1.1 

J034741.45+235818.7 

5154110 

78 


78 

56.9227179675443 

23.9718482388656 

2.6196128E-7 

2.6 

J034357.28+241320.4 

5146806 

94 


92 

55.9886710359267 

24.2223212907073 

1.8217466E-7 

1.8: 

J035041.20+253930.6 

5399682 

24 


24 

57.6716797131958 

25.6584928110123 

2.8959272E-7 

2.8' 

J034011.92+255232.3 

5374365 

60 


60 

55.0496801381824 

25.8756419920711 

2.8825687E-7 

2.8: 

J034412.14+235237.2 

5146567 

89 


89 

56.0506014764596 

23.8770056221031 

2.2743258E-7 

2.2' 

J035558.84+245740.0 

5294511 

60 


60 

58.9951560925242 

24.9611100221980 

1.9141787E-7 

1.9 

Table 4. DANCe / Tycho-2 catalogue of the Pleiades 

DANCe 

RA 

DEC 

pmRA 

pmDE e_pmRA e_pmDE BTmag e_BTmag VTmag 

e_V 

J033115.95+251519.7 

52.81646668 

25.25546943 

22.2 

-45.3 0.8 

0.8 8.44 

0.016 8.165 


J033128.34+214918.9 

52.86809418 

21.82190912 

23.8 

-48.8 0.7 

0.7 9.352 

0.021 8.884 


J033133.63+261555.8 

52.89010851 

26.26548691 

20.4 

-37.9 2.2 

2.4 11.62 

0.09 10.75 


J033305.82+220803.2 

53.27426485 

22.13422022 

21.4 

-49.3 1.6 

1.7 12.306 

0.225 11.521 


J033313.88+230023.1 

53.30782936 

23.00641262 

19.5 

-49.2 2.5 

2.6 13.42 

0.38 11.535 


J033401.81+245251.2 

53.50754459 

24.88089936 

20.8 

-42.1 1.6 

1.6 11.704 

0.079 10.932 


J033407.31+242040.0 

53.53045754 

24.34444364 

20.8 

-43.5 0.9 

0.9 10.119 

0.032 9.645 


J033458.66+233148.5 

53.74443215 

23.53014847 

20.2 

-42.5 1.0 

0.9 9.281 

0.022 8.891 


J033531.69+224924.9 

53.8820554 

22.82358994 

21.0 

-44.4 1.2 

1.3 10.484 

0.036 9.897 


J033608.46+272034.7 

54.03524149 

27.3429758 

22.5 

-48.2 1.3 

1.2 10.001 

0.028 9.367 


J033610.76+261552.4 

54.04483034 

26.26455153 

9.7 

-48.8 2.1 

2.1 10.942 

0.056 10.339 



